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The shortcomings of the Standard Big Bang Cosmological Model as well as their 
resolution in the context of inflationary cosmology are discussed. The inflation- 
ary scenario and the subsequent oscillation and decay of the inflaton field are 
then studied in some detail. The density perturbations produced during inflation 
and their evolution during the matter dominated era are presented. The tem- 
perature fluctuations of the cosmic background radiation are summarized. The 
non-supersymmetric as well as the supersymmetric hybrid inflationary model is 
introduced and the 'reheating' of the universe is analyzed in the context of the 
latter and a left-right symmetric gauge group. The scenario of baryogenesis via 
a primordial leptogenesis is considered in some detail. It is, finally, pointed out 
that, in the context of a supersymmetric model based on a left-right symmetric 
gauge group, hybrid inflation, baryogenesis via primordial leptogenesis and neu- 
trino oscillations are linked. This scheme, supplemented by a familiar ansatz for 
the neutrino Dirac masses and mixing of the two heaviest families and with the 
MSW resolution of the solar neutrino puzzle, implies that 1 eV ;S m„ T £ 9 eV. 
The mixing angle 6 MT is predicted to lie in a narrow range which will be partially 
tested by the Chorus/Nomad experiment. 



1 Shortcomings of the Big Bang Model 

The Standard Big Bang (SBB) Cosmological Model! has been very success- 
ful in explaining, among other things, the Hubble expansion of the universe, 
the existence of the Cosmic Background Radiation (CBR) and the abundances 
of the light elements which were formed during primordial nucleosynthesis. 
Despite its great successes, this model had a number of long-standing short- 
comings which we will now summarize: 

1.1 Horizon Problem 

The CBR, which we receive now, was emitted at the time of 'decoupling' of 
matter and radiation (which essentially coincides with the time of recombina- 
tion of atoms) when the cosmic temperature was w 3, 000 K. The decoupling 
time, td, can be calculated from 

To = 2.73 K = o(td) = ft_A 2/3 

T d 3,000 K a(t ) \t ) ' [ ' 

where to, Tq are the present cosmic time and temperature of CBR and a(t) 
is the dimensionless scale factor of the universe at cosmic time t normalized 
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so that a(io) = 1. It turns out that td ~ 200,000 h 1 years, where h is the 
present value of the Hubble parameter in units of 100 km sec -1 Mpc -1 . 

The distance over which the photons of the CBR have travelled since their 
emission is 

3< ~ 6,000 hr 1 Mpc , (2) 

which essentially coincides with the present particle horizon size. A sphere 
around us with radius equal to this distance is called the 'last scattering 
surface' since the CBR observed now has been emitted from it. The par- 
ticle horizon size at td was 2H~ 1 (td) = itd ~ 0.168 h^ 1 Mpc (H(t) being 
the Hubble parameter at cosmic time t) and expanded till now to become 
0.168 h- 1 (a(t )/a(t d )) Mpc w 184 h,- 1 Mpc. The angle subtended by this 
'decoupling' horizon at present is 9d ~ 184/6, 000 « 0.03 rads «2°. Thus, the 
sky splits into 4?r/(0.03) 2 « 14, 000 patches that never communicated causally 
before sending light to us. The question then arises how come the tempera- 
ture of the black body radiation from all these patches is so. accurately tuned 
as the measurements of the Cosmic Background Explorer □ (COBE) require 
(ST/Tk 6.6 x 10~ 6 ). 



f to dt' 



i- - 
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1.2 Flatness Problem 

The present energy density, p, of the universe has been observed to lie in the 
relatively narrow range 0.1p c $ p < 2p c , where p c is the critical energy density 
corresponding to a flat universe. The lower bound has been derived from 
estimates of galactic masses using the virial theorem whereas the upper bound 
from the volume expansion rate implied by the behavior of galactic number 
density at large distances. The Friedmann equation 

„ 8ttG k 8nG 
3 a z 3 

where H = a{t)/a{t) (overdots denote derivation with respect to cosmic time) 
is the Hubble parameter, G the Newton' s constant and k negative zero or pos- 
itive for an open, flat or closed universe respectively, implies that (/? — p c )/p c = 
3(8TrGp c )~ 1 (k/a 2 ) is proportional to a, for matter dominated universe. Con- 
sequently, in the early universe, |(p — p c )/pc\ ^ 1 and the question arises why 
the initial energy density of the universe was so finely tuned to be equal to its 
critical value. 
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1.3 Magnetic Monopole Problem 



This problem arises only if we combine the SBB Model with Grand Unified 
Theories!] (GUTs) of strong, weak and electromagnetic interactions. In accor- 
dance with GUTs, the universe underwent □ a phase transition during which 
the GUT gauge symmetry group, G, broke down to the Standard Model gauge 
group, Gs- This breaking was due to the fact that, at a critical temperature 
T c , an appropriate higgs field, 0, developed a non-zero vacuum expectation 
value (vev). Assuming that this phase transition was a second order one, we 
have (0)(T) « (4>){T = 0)(1 - T 2 /^ 2 ) 1 / 2 , m H (T) w A(0)(T), for the tempera- 
ture dependent vev and mass of the higgs field respectively at T < T c (A is an 
appropriate higgs coupling constant). 

The GUT phase transition produces magnetic monopoles El which are lo- 
calized deviations from the vacuum with radius ~ MZ , energy ~ Alx/ctc 
and = at their center (Mj is the GUT mass scale and ag = g G /4tt with 
gc being the GUT gauge coupling constant). The vev of the higgs field on a 
sphere, S* 2 , with radius 3> around the monopole lies on the vacuum man- 
ifold G/Gs and we, thus, obtain a mapping: S 2 — > G/Gs- If this mapping is 
homotopically non-trivial the topological stability of the magnetic monopole 
is guaranteed. 

Monopoles can be produced when the fluctuations of over = between 
the vacua at ±(0)(T) cease to be frequent. This happens when the free energy 
needed for to fluctuate from (0)(T) to zero in a region of radius equal to 
the higgs correlation length £(T) — mJ I 1 (T) exceeds T. This condition reads 
(47r/3)£ 3 AU & T, where AV ~ A 2 (0) 4 is the difference in free energy density 
between = and = (4>}{T). The Ginzburg temperature^, Tc, corresponds 
to the saturation of this inequality. So, at T < Tc, the fluctuations over = 
stop and (0) settles on the vacuum manifold G/Gs- At Tc, the universe splits 
into regions of size £g ~ (A 2 T C ) _1 , the higgs correlation length at Tc, with 
the higgs field more or less aligned in each region. Monopoles are produced at 
the corners where such regions meet (Kibble u mechanism) and their number 
density is estimated to be nyi ~ pCc? 3 ~ P^ 6 T{?, where p ~ 1/10 is a geo- 
metric factor. The 'relative' monopole number density then turns out to be 
tm = um/T ~ 10~ 6 . We can derive a lower bound on tm by pure causality 
considerations. The higgs field cannot be correlated at distances bigger than 
the particle horizon size, 2tc, at Tc- This gives the causality bound 

which implies that tm ~ 10~ 10 . 



3 



The subsequent evolution of monopoles, after Tq, is governed by the equa- 
tion El (overdots denote derivation with respect to cosmic time) 

driM 9 a 

= -Dn M - 3-tim , (5) 



dt 



where the first term in the right hand side (with D being an appropriate 
constant) describes the dilution of monopoles due to monopole-antimonopole 
annihilation while the second term corresponds to their dilution by the cosmo- 
logical expansion. The monopole-antimonopole annihilation proceeds as fol- 
lows. Monopoles diffuse towards antimonopoles through the plasma of charged 
particles, capture each other in Bohr orbits and eventually annihilate. The an- 
nihilation is effective provided the mean free path of the monopoles in the 
plasma of charged particles does not exceed their capture distance. This hap- 
pens at cosmic temperatures T > 10 12 GeV. The overall result is that, if 
the initial relative magnetic monopole density ru,m ~ 10~ 9 (;$ 10 -9 ), the final 
one rM,fin ~ 10 _9 (^ ?\M,in)- This combined with the causality bound yields 
fM,&n ~ 10 -10 . However, the requirement that monopoles do not dominate 
the energy density of the universe at nucleosynthesis gives 

r i\i(T m 1 MeV) < 10~ 19 , (6) 

and we have a discrepancy of about ten orders of magnitude. 

1.4 Density Fluctuations 

For structure formation i in the universe, we need a primordial density per- 
turbation, Sp/p, at all scales with a nearly flat spectrum 113. We also need 
some explanation of the temperature fluctuations, ST/T, of CBR observed by 
COBE0 at angles 9 Z 9d ~ 2 which violate causality (see Sec. 1.1). 
Let us expand 5p/p in plane waves 

S -P(x,t)=fd 3 kS- k (t)e^, (7) 
P J 

where a; is a comoving vector in 3-space and k is the comoving wave vector 
with k = \k\ being the comoving wave number (A = 2ir/k is the comoving wave 
length whereas the physical wave length is A p h ys = a(f)A). For A p h ys < H -1 , 
the time evolution of Sj, is described by the Newtonian equation 

4 + 2HS- k + V -^8- k = 4nG P S- k , (8) 
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where the second term in the left hand side comes from the cosmological ex- 
pansion and the third is the 'pressure' term (v s is the velocity of sound given 
by v 2 = dp/ dp, where p is the mean pressure). The right hand side of this 
equation corresponds to the gravitational attraction. 

For the moment, let us put H=0 (static universe). In this case, there 
exists a characteristic wave number fcj, the Jeans wave number, given by 
kj = AirGa 2 p/v 2 and having the following property. For k > kj, pressure 
dominates over gravitational attraction and the density perturbations just os- 
cillate, whereas, for k < fcj, gravitational attraction dominates and the density 
perturbations grow exponentially. In particular, for p—0 (cold dark matter), 
v s = and all scales are Jeans unstable with 

5 U cx exp(t/r) , r = (47rGp)~ 1/2 . (9) 

Now let us take H ^ 0. Since the cosmological expansion pulls the particles 
apart, we get a smaller growth: 

5- k oc a{t) oc t 2/3 , (10) 

in the matter dominated case. For a radiation dominated universe {p =^ 0), 
we get essentially no growth of the density perturbations. This means that, in 
order to have structure formation in the universe, which requires Sp/ p ~ 1, we 
must have 

(^)e,u~4x l(T 5 (^r 2 , (11) 

at the equidensity point (where the energy densities of matter and radia- 
tion coincide), since the available growth factor for perturbations is given by 
tto/ciequ ~ 2.5 x lO A (Q a h) 2 . Here ilo = po/p c , where po is the present energy 
density of the universe. The question then is where these primordial density 
fluctuations originate from. 

2 Inflation 

InflationEU'El is an idea which solves simultaneously all four cosmological puz- 
zles and can be summarized as follows. Suppose there is a real scalar field (f> 
(the inflaton) with (symmetric) potential energy density V((f>) which is quite 
'flat' near = and has minima at cf> = ±(</>) with V(±(<fi)) = 0. At high 
enough T" s, <j) — in the universe due to the temperature corrections in V{<p). 
As T drops, the effective potential density approaches the T=0 potential but 
a little potential barrier separating the local minimum at <j) — and the vacua 
at (f> = ±{4>) still remains. At some point, <f> tunnels out to 4>\ <C ((f)) and a 
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bubble with = <pi is created in the universe. The field then rolls over to the 
minimum of V(<p) very slowly (due to the flatness of the potential). During 
this slow roll-over, the energy density p sa V(<t> — 0) = Vb remains essentially 
constant for quite some time. The Lagrangian density 



L = - V{4>) (12) 



gives the energy-momentum tensor 

T^ v = -3^8^ + 8^ {^dx^^-V^ , (13) 

which during the slow roll-over takes the form v w — Vq 5^ v . This means 
that p w — p w Vo, i.e., the pressure p is negative and equal in magnitude 
with the energy density p, which is consistent with the continuity equation 
p = —3H(p + p). Since, as we will see, a(t) grows very fast, the 'curvature' 
term, fc/a 2 , in Eqj^ becomes subdominant and we get 

H^( a -X= 8 -^V , (14) 



which gives a(t) oc e Ht , H 2 = (8ttG/3)Vq — constant. So the bubble expands 
exponentially for some time and a(t) grows by a factor 

^4 = e*pH(t f - U) = expiJr , (15) 

(X[ti ) 

between an initial (ti) and a final (tf) cosmic time. _ 

The inflationary scheme just described, which is known as the newt3 in- 
flationary scenario (with the inflaton field starting from the origin, 0=0), is 
certainly not the only realization of the idea of inflation. Another interest- 
ing possibility is to consider the universe as it emerges at the Planck time 
t P = mp l ,m P = Mp/VSn with M P = 1.22 x 10 19 GeV, where the fluctu- 
ations of gravity cease to exist. We can imagine a region of size lp w nip 1 
where the inflaton field acquires a large and almost uniform value and carries 
negligible kinetic energy. Under certain circumstances this region can inflate 
(exponentially expand) as <p rolls down towards its vacuum value. This type 
inflation with the inflaton starting from large values is known as the chaoticE 
inflationary scenario. 

We will now show that, with an adequate number of e-foldings, N = 
Ht, the first three cosmological puzzles are easily resolved (we will leave the 
question of density perturbations for later). 
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2.1 Resolution of the Horizon Problem 

The particle horizon during inflation (exponential expansion) 

r f dt 1 

d W = eHt "W ~ H^e^Hit - U) , (16) 
Jti e 

for t — ti^$> H~ x , grows as fast as a{t). At the end of inflation (t = tf), d(tf) rs 
H~ 1 ex.pHr and the field <fi starts oscillating about the minimun of the potential 
at cj) = (</>). It then decays and 'reheats'Ej the universe at a temperature T r ~ 
10 9 GeV (see Sec.[lO|). The universe, after that, goes back to normal big bang 
cosmology. The horizon d(tf) is stretched during the period of ^-oscillations 
by some factor ~ 10 9 depending on details and between T r and the present 
era by a factor T r /T . So it finally becomes equal to H - 1 e Hr 10 9 (T r /T ), 
which should exceed 27?^ 1 in order to solve the horizon problem. Taking 
V w M%, M x ~ 10 16 GeV, we see that, with N — Ht > 55, the horizon 
problem is evaded. 

2.2 Resolution of the Flatness Problem 

The curvature term of the Friedmann equation, at present, is given by 

k fk\ ?Ht m flQ- 13 GeV\ 2 . . 

where the terms in the right hand side correspond to the curvature term before 
inflation, and its growth factors during inflation, during <fi -oscillations and after 
'reheating' respectively. Assuming (k/a 2 )i,i ~ (8irG/3)p ~ H 2 (p Vq), we 
get k/alH 2 - 10 48 e~ 2IiT which gives (p - p c )/p c = il - 1 = k/alHg < 1, 
for Ht 3> 55. In fact, strong inflation implies that the present universe is flat 
with a great accuracy. 

2.3 Resolution of the Monopole Problem 

It is obvious that, with a number of e- foldings £ 55, the primordial monopole 
density is diluted by at least 70 orders of magnitude and they become totally 
irrelevant. Also, since T r <C tum, there is no production of magnetic monopoles 
after 'reheating'. 
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3 Detailed Analysis of Inflation 

The Hubble parameter is not exactly constant during inflation as we, naively, 
assumed so far. It actually depends on the value of (f>: 

H\4>) = . (18) 

To find the evolution equation for (f> during inflation, we vary the action 

S = J v/^dct^) d A x (-d^t - V(cf>) + M{4>Yj , (19) 

where g is the metric tensor and M{<f)) represents the coupling of <fr to 'light' 
matter causing its decay. We find 

^ + 3H4> + r^ + v'{4) = o , (20) 

where-the prime denotes derivation with respect to <p and is the decay 
width Ii3 of the inflaton. Assume, for the moment, that the decay time of cf>, 
td = rT , is much greater than if -1 , the expansion time for inflation. Then 

the term T^cj) can be ignored and Eqj2^ reduces to 

4> + 3H<p + V'((j)) = . (21) 

Inflation is by definition the situation where <f> is subdominant to the 'friction' 
term ?>H <\> in this equation (and the kinetic energy density is subdominant to 
the potential energy density). Eq.|^ then further reduces to the inflationary 
equation!!!! 

3H(j>= -V'(^) , (22) 



which gives 

~ v"(<t>)j> , V'(4>) 



ZH{4>) ZH 2 {4>) 



H'(<f>)j> . (23) 



Comparing the two terms in the right hand side of this equation with the 
'friction' term in Eq.^l], we obtain the conditions for inflation (slow roll condi- 
tions): 

V"(<t>) 



_ IVlp 

8tt 



V{4>) 



The end of the slow 'roll-over' occurs when either of the these inequalities is 
saturated. If 0/ is the value of 4> at the end of inflation, then tf ~ H~ x (</>/). 
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The number of e-foldings during inflation can be calculated as follows: 



, s „ fa(t f )\ f tf TT1 f+f H[6) 71 t*' 2,H 2 U)d(j) 

(25) 

where Eqs.|f5|, ^2| and the definition of H = a/a were used. For simplicity, we 
can shift the field <fi so that the global minimum of the potential is displaced 
at (j> = 0. Then, if V(<fi) = X(f) 1 ' during inflation, we have 

+' 3tf 2 (0)# Q „ [+' V(<f>)d<f> 4nG i 

(26) 



Assuming that <pi 3> (j>f, this reduces to N(<f>) = (AnG/v)(j) 2 . 
4 Coherent Field Oscillations 

After the end of inflation at cosmic time tf, the term <f> takes over and Eq.pT] 
reduces to <j> + V'(<p) — 0, which means that <f> starts oscillating coherently 
about the global minimum of the potential. In reality, due to the 'friction' 
term, <fi performs damped oscillations with a rate of energy density loss given 
by 

d 4 = i ( X J 2 + V {</>)) = SH^ = -3H(p + p) , (27) 



dt dt\2 

where p = <f> 2 jl + V(<p) and the pressure p — cb 2 /2 — V ((/)). Averaging p over 
one oscillation of <f>, we writeES p + p = "/p. Eq.^7| then becomes p = —3Hjp, 
which gives dp/ p = —3jda/a and p cx a~ 37 . The Friedmann equation then 
takes the form a/a oc a~ 37 ' 2 and we obtain a(t) oc i 2 ' 37 . 

The number 7 for an oscillating field can be written as (assuming a sym- 
metric potential) 

fipdt st^ip/m 



7= £ 7t— = 7T7T, ■ (28) 



where T and max are the period and the amplitude of the oscillation respec- 
tively. From the equation p — <f> 2 /2 + V(<fi) — V m ^, where V max is the maximal 
potential energy-density, we obtain <fi = y/2(V max — V ((/>)). Substituting this 
in Eq.^H we getliM 

= 2/^(1-7^)^ 
J^— (l_y/y max) -i/2 # 
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For a potential of the simple form V(4>) — Xtp 1 " , 7 is readily found to be given by 
7 = 2v/{v+2). Consequently, in this case, p oc a ~ 6,y/(,y+2) and a(t) oc t^ +2 ^ 3l/ . 
For v = 2, in particular, one has 7=1, p oc a~ 3 , a(t) oc f 2 / 3 and the oscillating 
field behaves like pressureless 'matter'. This is not unexpected since a coherent 
oscillating massive free field corresponds to a distribution of static massive 
particles. For ^=4, however, we obtain 7 = 4/3, p oc a -4 , a(t) oc t 1 / 2 and the 
system resembles 'radiation'. For v = 6, one has 7 = 3/2, p oc a -4 5 , a(t) oc 
i 4 / 9 and the expansion is slower than in a 'radiation' dominated universe (the 
pressure is higher than in 'radiation'). 

5 Decay of the Field </> 



Reintroducing the 'decay' term Y^cj), Eq,20 can be written as 



dt V 2 

which is solved Bill by 



p = i(U 2 + v W) =-(3# + r )0 2 , (801 



pit) = Pf (^) 37 exp[- 7 rv; - t f )] , (31) 

where pf is the energy density at the end of inflation at cosmic time tf. The 
second and third factors in the right hand side of this equation represent the 
dilution of the field energy due to the expansion of the universe and the decay 
of <fr to light particles respectively. 

All pre-existing 'radiation' (known as 'old radiation') was diluted by infla- 



tion, so the only 'radiation' present is the one produced 
is known as 'new radiation'. Its energy density satisfies 



L5LP-I 



ie decay of <f> and 
the equation 



p r = —AHp r + 7 r p , (32) 

where the first term in the right hand side represents the dilution of radiation 
due to the cosmological expansion while the second one is the enersx-xLensity 
transfer from tfi to 'radiation'. Taking p r (tf)=0, this equation givesEj£3 

For tf <^td and v = 2, this expression is approximated by 

/ . \ -8/3 At / ,1 \ 2/3 

<* w ="(£) 1(h) e ' T ' l ' dt '- (34) 
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which, using the formula 



dx = er u 



p+k 



fc=0 



p(p + l)---(p + fe) ' 



can be written as 



Pr = - P IV 
5 



i + |iV + A(iy) a 



(35) 



(36) 



with p = pf(t/tf)~ 2 exp(— F$t) being the energy density of the field <j> which 
performs damped oscillations and decays into 'light' particles. 

The energy density of the 'new radiation' grows relative to the energy 
density of the oscillating field and becomes essentially equal to it at a cosmic 
time td — TT 1 as one can deduce from Eq.36. After this time, the universe 



enters into the radiation dominated era and the normal big bang cosmology is 
recovered. The temperature at td, T r (td), is historically called the 'reheating' 
temperature although no supercooling and subsequent reheating of the universe 
actually takes place. Using the time to temperature relation for a radiation 
dominated universe we find that 



45 



\16 



1/4 



(37) 



where is the effective number of degrees of freedom. For a potential of the 
type V((p) — \<t> v , the total expansion of the universe during the period of 
damped field oscillations is 



a(td) 
«(*/) 



(38) 



6 Density Perturbations 

We are now ready to sketch how inflation solves the density fluctuation problem 
described in Sec. 



1.4 



As a matter of fact, inflation not only homogenizes the 
universe but also provides us with the primordial density fluctuations necessary 
for the structure formation in the universe. To understand the origin of these 
fluctuations, we must first introduce the notion of 'event horizon'. Our 'event 
horizon', at a cosmic time t, includes all points with which we will eventually 
communicate sending signals at t. The instantaneous (at cosmic time t) radius 
of the 'event horizon' is 

f°° dt' 

d e (t)=a(t) J ^ ■ (39) 
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It is obvious, from this formula, that the 'event horizon' is infinite for matter 
or radiation dominated universe. For inflation, however, we obtain a slowly 
varying event horizon with radius d e (t) = H^ 1 < oo. Points, in our 'event 
horizon' at t, with which we can communicate sending signals at i, are even- 
tually pulled away by the 'exponential' expansion and we cease to be able to 
communicate with them again emitting signals at later times. We say that 
these points (and the corresponding scales) crossed outside the event horizon. 
The situation is very similar to that of a black hole. Indeed, the exponentially 
expanding (de Sitter) space is like a black hole turned inside out. This means 
that we are inside and the black hole surrounds us from all sides. Then, exactly 
as in a black hole, there are quantun^fiuctuations of the 'thermal t#H^ gov- 
erned by the 'Hawking temperature' li*3'E3 Th = H/2n. It turns outEirE3 that 
the quantum fluctuations of all massless fields (the inflaton is nearly massless 
due to the 'flatness' of the potential) are 5(fi = H/2tt = Th- These fluctua- 
tions of <f> lead to energy density fluctuations 8p = V'((f))6<fi. As the scale of 
this perturbations crosses outside the event horizon, they become Ea classical 
metric perturbations. 

The evolution of these fluctuations outside the 'inflationary horizon' is 
quite subtle and involved due to the gauge freedomJii general relativity. How- 
ever, there is a simple gauge invariant quantity Ej C w Sp/(p + p), which 
remains constant outside the horizon. Thus, the density fluctuation at any 
present physical (comoving) scale £, (5p/p)e, when this scale crosses inside the 
'post-inflationary' particle horizon (j>—0 at this instance) can be related to the 
value of £ when the same scale crossed outside the inflationary event horizon 
(symbolically at £ ~ H^ 1 ). This latter value of C can be found using Eq.p2| 
and turns out to be 



Taking into account an extra 2/5 factor from the fact that the universe is 
matter dominated when the scale £ re-enters the horizon, we obtain 



The calculation of (fig, the value of the inflaton field when the 'comov- 
ing' scale £ crossed outside the event horizon, goes as follows. A 'comov- 
ing' (present physical) scale £, at T r , was equal to £(a(td) / a(toj) = £(T /T r ). 
Its magnitude at the end of inflation (t = tf) was £(To/T r )(a(tf)/a(t ( i)) = 
£{T /T r )(t f /t d )^ +2 ^' iu = £ phys {t f ), where the potential V(cfi) = X(fi v was as- 
sumed. The scale £, when it crossed outside the inflationary horizon, was equal 




(40) 




16^6^ V 3 / 2 (<fi e ) 
5 M%V'(<t> t ) 



(41) 
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to H 1 (<fre). We, thus, obtain 

H-\^)e N ^ =i phys (t f ) . (42) 

Solving this equation, one can calculate <pi and, subsequently, N(<fii) = Ni, 
the number of e-foldings the scale I suffered during inflation. In particular, for 
our present horizon scale t s» 2iJ ~ 1 ~ 10 4 Mpc (Hq is the present value of the 
Hubble parameter), it turns out that Nh ~ 50 — 60. 



Now, taking the potential V(<fi) — A0 4 , Eqs.|26j, [ll| and |42| give 



'V\ - ! \ ( ' 7r ^i/2 ( <k V 4^ A i/ 2 (NA 3/2 _ (43) 



p / « 5 \Mp / 5 

The measurements of COBe!, {Sp/p)H a ~ 6 x 10~ 5 , then imply that A w 
6 x 10" 14 for Nh ~ 55. Thus, we see that the inflaton must be a very weakly 
coupled field. In non-supersymmetric GUTs, the inflaton must necessarily be 
a gauge singlet field since otherwise radiative corrections will certainly make it 
strongly coupled. This is, undoubtedly, not a very satisfactory situation since 
we are forced to introduce an otherwise unmotivated extra ad hoc very weakly 
coupled gauge-singlet. In supersymmetric GUTs, however, the inflaton could 
be identified c3 with a conjugate pair of gauge non-singlet fields 0, <j>, already 
existing in the theory and causing the gauge symmetry breaking. Absence of 
strong radiative corrections from gauge interactions is guaranteed, in this case, 
by the mutual cancellation of the D-terms of these fields. 

The spectrum of the density fluctuations emerging from inflation can also 
be analyzed. For definiteness, we will again take the potential V(4>) = Xcjf . 
One then finds that (5p/p)i is proportional to <j)^ /+2 ^ 2 which, combined with 



the fact that N{4>g) is proportional to 4> 2 (see Eq.26), gives 

Sp\ I Sp\ ( Ni 



P J l V P J H \ N Ho 



(44) 



The scale I divided by the size of our present horizon (« 10 4 Mpc) should equal 
exp(JVi - N Ha ). This gives N e /N Ha = 1 + fc^/lO 4 ) 1 /^ which expanded 
around I rj 10 4 Mpc and substituted in Eq.[l4| yields 



6p\ ( bp 



PJi \P J ffr ,Vl0 4 Mpc 



(45) 



with a s = (y + 2)/4iV# . For v = 4, a s w 0.03 and the fluctuations are 
essentially scale independent. 
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7 Density Fluctuations in 'Matter' 

We will now discuss the evolution of the primordial density fluctuations after 
their scale enters the post-inflationary horizon. To this end, we introduce 
the 'conformaP time, 77, so that the Robertson- Walker metric takes the form 
of a conformally expanding Minkowski space: 

ds 2 = -dt 2 + a 2 {t) df 2 = a 2 (ri) {-drj 2 + df 2 ) , (46) 

where f is a comoving 3-vector. The Hubble parameter now takes the form 
H = a(t)/a(t) — a' (rj) / a 2 (rj) and the Friedmann equation can be rewritten as 

1 f a'\ 2 8ttG 

= , (47) 



,2 



a 



where primes denote derivation with respect to the 'conformal' time rj. The 
continuity equation takes the form p' = —3H(fi + p) with H = a' /a. For a 
matter dominated universe, p oc a~ 3 which gives a = (77/770) 2 an d a' /a = 2/rj 
(r]o is the present value of rj). 

The 'Newtonian' Eq.|§| can now be written in the form 

S 'iM + ^S'kiv) - ^Gpa 2 S- k {r,) = , (48) 

and the grxaving (Jeans unstable) mode 5^ (77) is proportional to ir 2 and can be 
expressed E3 as 

' krj^ 2 



h(v) = eH [i:) J ( fc ) . ( 49 ) 

where s(k) is a Gaussian random variable satisfying 

<s(fc)>=0, < s(k)s(k') >= ^5(k - k') , (50) 

k 

and en is the amplitude of the perturbation when its scale crosses inside the 
'post-inflationary horizon' . The latter can be seen as follows. A 'comoving' 
(present physical) length I crosses inside the post-inflationary horizon when 
a£/2n = if" 1 = a 2 /a' which gives £/2vr = fc" 1 = a/ a' = r) H /2 or kr) H /2 = 1, 
where ijh is the conformal time at horizon crossing. This means that, at 
horizon crossing, 8 k {rjH) = ^Hs(k). For scale invariant perturbations, the 
amplitude en is constant. The gauge invariant perturbations of the 'scalar 
gravitational potential' are giveno by the Poisson's equation, 

$ = -Ai,G^p8- k (rf) . (51) 
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From the Friedmann Eq.|47j, we then obtain 

$ = -^Hs{k) • (52) 

The spectrum of the density perturbations can be characterized by the 
correlation function 

Z(f)=<6*(x,r,)6(x + f,ri)> , (53) 

where 

~6(x )V ) = Jd 3 kS- k ( v )e^ . (54) 
Substituting Eq.^9| in Eq.^3] and using Eq.|5(i|, we obtain 

«f) = / rf3fce ^(y) 4 ^' (55) 

and, thus, the spectral function P(k,rj) = e 2 H (rf /\&)k is proportional to k for 
eh constant. We say thaL-an this case, the 'spectral index' n = 1 and we 
have a Harrison- Zeldovicht2l flat spectrum. In the general case, P cx k n with 
n = l- 2a s (see Eq.||). For V{<j>) = A0 4 , we get n w 0.94. 

8 Temperature Fluctuations 

The density inhomogeneities produce temperature fluctuations in the CBR. 
For angles 9 > 2 °, the dominant effect is the scalar Sachs- Wolfe cB effect. 
Density perturbations on the 'last scattering surface' cause scalar gravitational 
potential fluctuations, <j>, which, in turn, produce temperature fluctuations 
in the CBR. The physical reason is that regions with a deep gravitational 
potential will cause the photons to lose energy as they climb up the well and, 
thus, appear cooler. For 8^2°, the dominant effects are: i) Motion of the 
last scattering surface causing Doppler shifts, and ii) Intrinsic fluctuations of 
the photon temperature, T 7 , which are more difficult to calculate since they 
depend on microphysics, the ionization history, photon streaming and other 
effects. 

The temperature fluctuations at an angle 6 due to the scalar Sachs- Wolfe 
effect turn outE3 to be given by (ST/T)g = — $^/3, where t is the 'comov- 
ing' scale on the 'last scattering surface' which subtends the angle 9 [ £ w 
100 /i _1 (6'/degrees) Mpc ] and $£ the corresponding scalar gravitational po- 
tential fluctuations. From Eq|3^, we then obtain (6T/T)$ = (eh /2)s(k), which 



15 



using Eq.^ becomes 

(?) 



(56) 



The COBE scale (our present horizon )corresponds to 
and give 



60 °. EqsM 



41 



ST 



dp 



P J t 



M%V'{<t>t) 



oc n: 



(57) 



Analyzing the temperature fluctuations in spherical harmonics we can obtain 
the quadrupole anisotropy due to the scalar Sachs- Wolfe effect, 



ST 

1 / Q-S 

For V(4>) — \(f)" , this becomes 



/32^\ 1/2 V 3 ' 2 {^) 



V 45 y 



l) 



(58) 



ST 



32tt 
~45~ 



1/2 



32tt\ 1/2 A 1 / 2 /i/A^ 



45 



An 



n: 



(59) 

Comparing this with the COBE0 measurements, (5T/T)q « 6.6 x 10 6 , we 
obtain A ~ 6 x 10~ 14 , for v = 4 and number of e- foldings suffered by our 
present horizon scale during inflation N^ H -i = Nq w 55. 

There are also 'tensor' E3 (gravitational wave) fluctuations in the temper- 
ature of CBR. The quadrupole tensor anisotropy is 



ST\ 
~T) 



0.77 



Q-T 



M 2 



The total quadrupole anisotropy is given by 



ST 



ST 
T 



and the ratio 



(ST/T) 2 Q _ T 
(ST/T)%_ S 



)-s 



0.27 



ST 
T 



-1 1/2 



M P V>(& 
V(<j>e) 



(60) 



(61) 



(62) 



For V{<f>) — \(f> v , we obtain r w 3.4 v/Nh *C 1, and the 'tensor' contribution 
to the temperature fluctuations of the CBR is negligible. 
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9 Hybrid Inflation 



9.1 The non Super 'symmetric Version 

The most important disadvantage of the inflationary scenarios described so far 
is that they neecLextremely small coupling constants in order to reproduce the 
esults of COBE0. This difficulty was overcome some years ago by LindeEll who 
proposed, in the context of non-supersymmetric GUTs, a clever inflationary 
scenario known as hybrid inflation. The idea was to use two real scalar fields 
X and a instead of one that was normally used. The field x provides the 
vacuum energy which drives inflation while a is the slowly varying field during 
inflation. The main advantage of this scenario is that it can reproduce the 
observed temperature fluctuations of the CBR with 'natural' values of—the 
parameters in contrast to previous realizations of inflation (like the 'new' 113 or 
'chaotic' li-3 inflationary scenarios). The potential utilized by Linde is 

V( X ,«) = k 2 (m>-4)\^ + ^ , (63) 



where k, A are dimensionless positive coupling constants and M, m mass 
parameters. The vacua lie at (x) — ±2M, (a) = 0. Putting m=0, for the 
moment, we observe that the potential possesses an exactly flat direction at 
X = with V(x = 0, er) = k 2 M 4 . The mass squared of the field x along this 
flat direction is given by m 2 , = —k 2 M 2 + \\ 2 u 2 and remains non-negative for 
cr > cr c = V2kM/X. This means that, at x = and a > cr c , we obtain a valley 
of minima with flat bottom. Reintroducing the mass parameter m in Eq.p3|, 
we observe that this valley acquires a non-zero slope. A region of the universe, 
where x and a happen to be almost uniform with negligible kinetic energies and 
with values close to the bottom of the valley of minima, follows this valley in 
its subsequent evolution and undergoes inflation. The quadrupole anisotropy 



of CBR produced during this inflation can be estimated, from Eq.58, to be 
'8T\ /16tt\ 1/2 Xk 2 M 5 



- 



T J Q ~ V 45 J Mlm 2 ' (64) 

The COBE I result, (5T/T) Q » 6.6 x 1(T 6 , can then be reproduced with 
M w 2.86 x 10 16 GeV (the supersymmetric GUT vev) and m « 1.3 kVA x 10 15 
GeV ~ 10 12 GeV for k, A ~ 10~ 2 . Inflation terminates abruptly at a = <j c 
and is followed by a 'waterfall', i.e., a sudden entrance into an oscillatory 
phase about a global minimum. Since the system can fall into either of the 
two available global minima with equal probability, topological defects are 
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copiously produced if they are predicted by the particular particle physics 
model one is considering. 

9.2 The Supersymmetric Version 

The hybrid inflationary scenario is 'tailor made' for application to super- 
symmetric GUTs except that the mass of a, m, is unacceptably large for su- 
persymmetry, where all scalar fields acquire masses of order ms ~ 1 TeV from 
soft supersymmetry breaking. To see this, consider a supersymmetric GUT 
with a (semi-simple) gauge group G of rank > 5 with G — > Gs (the Standard 
Model gauge group) at a scale M ~ 10 16 GeV. The spectrum of the theory 
below M is assumed to coincide with the Minimal Supersymmetric Standard 
Model (MSSM) spectrum plus Standard Model singlets so that the successful 
predictions for a s , sm 2 6w are retained. The theory may also possess global 
symmetries. The breaking of G is achieved through the superpotential 

W = kS(-M 2 + #), (65) 

where <f>, <j> is a conjugate pair of Standard Model singlet left handed superfields 
which belong to non-trivial representations of G and reduce its rank by their 
vevs and S is a gauge singlet left handed superfield. The coupling constant K 
and the mass parameter M can be made real and positive by suitable redefini- 
tions of the phases of the superfields. This superpotential is the most general 
renormalizable superpotential consistent with a U(l) R symmetry under which 
W -> e l6 W, S -> e i9 S, 4><j> -> # and gives the potential 

V = K 2 \M 2 -^\ 2 +K 2 \S\ 2 (\<j>\ 2 + \4>\ 2 ) 

+ D — terms. (66) 

Restricting ourselves to the D-flat direction <fi* = (f> which contains the super- 
symmetric minima and performing appropriate gauge and R transformations, 
we can bring S, <f>, <f) on the real axis, i.e., S = <f> = <j) = x/2, where 

cr, x are normalized real scalar fields. The potential then takes the form in 
Eq.|63| with k = A and m = and, thus, Linde's potential for hybrid inflation 
is almost obtainable from supersymmetric GUTs but without the mass term 
of a which is, however, of crucial importance since it provides the slope of the 
valley of minima necessary for inflation. 

One way to obtain a valley of minima useful for inflation iso to replace 
the renormalizable trilinear term in the superpotential W in Eq.^5| by the next 
order non-renormalizable coupling. Another way, which we will adopt here, 
is E3 to keep the renormalizable superpotential in Eq.pq and use the radiative 
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corrections along the inflationary valley (cf> — cj) = 0,S>S c = M). In fact, 
due to supersymmetry breaking by the 'vacuum' energy density , k 2 M 4 , along 



this valley, there are important radiative corrections. At one 
sufficiently larger than S c , the inflationary potential is given t 



V eS (S) = k 2 M 4 



1 + t in 



16tt 2 V V A 2 / 2 12S 4 



3 



ap, and for S 
by 



(67) 



where A is a suitable mass renormalization scale. Using this effective poten- 
tial and Eq.|5^, one finds that the cosmic microwave quadrupole anisotropy 
(ST/T) Q » 8n(N Q /A5) 1 / 2 (x Q /y Q )(M/M P ) 2 . Here y Q = x Q (l - 7/(12x 2 ) + 
■ ■ •) with xq = Sq /M, and Sq is the value of the scalar field S when the scale 
which evolved to the present horizon size crossed outside the de Sitter (infla- 
tionary) horizon. Also from Eq.^37|, one finds k m (87r 3 / 2 /y]VQ) y Q (M/Mp). 

The inflationary phase ends as S approaches S c from above. Writing S = 
xS c , x = 1 corresponds to the phase transition from G to G$ which, as it 
turns out, more or less coincides with the end of the inflationary phase (this is 
checked by noting the amplitude of the quantities e and r\ in Eq.p4[). Indeed, the 
50 — 60 e- foldings needed for the inflationary scenario can be realized even with 
small values of xq. For definiteness, we will take xq w 2 (see SecJT^). From 
COBE0 one then obtains M w 5.5 x 10 15 GeV and k«4.5x 10" 3 for Nq m 56. 
Moreover, the primordial density fluctuation spectral index n ~ 0.98. We see 
that the relevant part-piipflation takes place at S ~ 10 16 GeV. An important 
consequence of this isoaa that the supergravity corrections can be negligible. 

In conclusion, it is important to note that the superpotential W in Eq.|65| 
leads to hybrid inflation in a 'natural' way. This means that a) there is no 
need of very small coupling constants, b) W is the most general renormalizable 
superpotential allowed by the gauge and R symmetries, and c) supersymmetry 
guarantees that the radiative corrections do not invalidate inflation. They 
rather provide a slope along the inflationary trajectory which drives the inflaton 
towards the supersymmetric vacua. 



10 'Reheating' in Supersymmetric Hybrid Inflation 

In order to discuss the 'reheating' of the universe after inflation, we need to 
be a little more specific about the underlying particle physics model. To this 
end, we considertj a supersymmetric model based on the left-right symmetric 
gauge group G LR = SU(3) C x SU(2) L x SU(2) R x U(1) B -l- Of course, it is 
anticipated that Glr is embedded in a grand unified theory such as 5*0(10) 
or SU(3) C x SU(3) L x SU(S)r. The breaking of Glr to G s is achieved by the 
renormalizable superpotential in Eqj6q, where 0, 4> are now identified with the 
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Standard Model singlet components of a conjugate pair of SU(2)r x U(1)b-l 
doublet left handed superfields. 

The inflaton (oscillating system) consists of the two complex scalar fields 
S and 9 — (8<j> + 8<j>)/y/2, where 8(f) = <j> — M, 8(f> — — M, with mass 
Triinfi =\[2kM. We mainly concentrate on the decay of 9. Its relevant coupling 
to 'matter' is provided c3 by the non-renormalizable superpotential coupling 
(in symbolic form) 

1 (M v c 



2 V M 2 < < 68 > 
where M v o denotes the Majorana mass of the relevant right handed neutrino 
v c . Without loss of generality, we assume that the Majorana mass matrix of 
the right handed neutrinos has been brought to diagonal form with positive en- 
tries. Clearly, 9 decays predominantly into the heaviest right handed neutrino 
permitted by phase space, i.e., with M„c < mj„/z/2. _ . 

The field S is not important for reheating since it can rapidly decayEJ into 
higgsinos through the renormalizable superpotential term ^Sh^'h^ allowed 
by gauge symmetry, where denote the electroweak higgs doublets 

which couple to the up and down type quarks respectively, and £ is a suit 
coupling constant. Note that, after supersymmetry breaking, S acquires! 
an expectation value (S) ~ ms (ms ~ 1 TeV being the magnitude of super- 
symmetry breaking in the visible sector) and generates the \x term. 

Following standard procedures (see Eq. |3?]), and assuming the MSSM spec- 
trum, the 'reheating' temperature T r is found to be given by 

T r w i(IW P ) 1/2 , (69) 

where T g w (1/16tt)(\/2A/„c /M) 2 V2kMJs the decay rate of 9. Substituting k 
as a function of Nq, uq and M, we finda 



1 / 56 N 1/4 



12 \N, 



fy^M v *. (70) 



For iq»2 and Nq « 56, we have T*-£i M u o /9.23. We will assume that T r is 
restricted by the gravitino constraint c3, T r < 10 9 GeV. Note that T r is closely 
linked to the mass of the heaviest v c satisfying M„c < min fi/l. 



11 Baryogenesis via Leptogenesis 

11.1 Primordial Leptogenesis 

In the hybrid inflationary models under consideration here, it is not conve- 
nient to produce the baryon asymmetry of the universe (BAU) in the usual 
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way, i.e., through the decay of color 3, 3 fields (g, g c ). Some of the reasons 
are the following: i) For theories where leptons and quarks belong to different 
representations of the unifying gauge group G (which is the case, for exam- 
ple, for G = Glr or SU(3) C x SU(3)l x SU(3)r), the baryon number can 
be made almost exactly conserved by imposing an appropriate discrete sym- 
metry. In particular, for G — Glr, we can impose E3 a discrete symmetry 
under which q — > —q, q c — > — q c , q — > —q, q c — > — q c and all other superfields 
remain invariant (q, q c , q, q c are superfields with the quantum numbers of the 
quarks, antiquarks and their conjugates respectively), ii) For theories where 
such a discrete symmetry is absent, we could, in principle, use as inflaton a 
pair of conjugate Standard Model singlet superfields N, N which decay into 
g, g c . For G = SU(3) C x SU{3) L x SU(3) R , for example, N (N) could be 
the Standard Model singlet component of the (1, 3, 3) ( (1, 3, 3) ) super- 
fields with zero U(1)b-l charge. But this is again unacceptable since the 
breaking of SU(3) C x SU(3) L x SU(3) R by the vevs of N, N predicts^ mag- 
netic monopoles which can then be copiously produced after inflation. Also 
T r < 10 9 GeV (gravitino constraint E3) implies m g < 10 10 GeV ( from the 
coupling (m g / (N))Ngg c ) which leads to strong deviation from MSSM and 
possible proton decay problems. _ 

So it is preferable to produce first a primordial lepton asymmetry l3 which 
can then be partially turned into the observed baxyon asymmetry of the uni- 
verse by the non-perturbative 'sphaleron' effects El of the electroweak sector. 
In the particular model based on Glr which we consider here, this is the 
only way to produce the BAU since the inflaton decays into right handed neu- 
trinos. Their subsequent decay into ordinary higgs particles (higgsinos) and 
light leptons (sleptons) can produce the primordial lepton asymmetry. It is 
important, though, to ensure that this primordial lepton asymmetry is not 
erased by lepton number violating 2 — ► 2 scatterings such as 11 — > * 
or lh^ — > Ths 1 ' * (I represents a lepton doublet) at all temperatures between 
T r and lQfLGeV. This is automatically satisfied since the lepton asymmetry is 
protectedcj by supersymmetry at temperatures between T r and T ~ 10 7 GeV, 
and for T < 10 7 GeV, these 2 —> 2 scatterings areCJ well out of equilibrium. 

The lepton asymmetry produced by the out-of-equilibrium decay (M„c ^> 
T r ) of the right handed neutrinos vf, which emerged from the inflaton decay, 
is calculated to beEJeJ 

S~ 167rmmfl 2^^ rh >\ {h (l))\2 {U M D, M D,i > W 

where til and s are the lepton number and entropy densities, M D 1 is the 
diagonal 'Dirac' mass matrix, U a unitary transformation so that UM D ' is 
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the 'Dirac' mass matrix in the basis where the 'Majorana' mass matrix of v c 's 
is diagonal (see Sec.p"2|) and | {/t.^ 1 ^ } | 174 GeV for large tan/3. The function 



with g(r) ~ 1/r as r — > oo. Here we took into account the following prefactors: 
i) At 'reheating', we have ni n fimi n fi = (7r 2 /30)g*T^ (ni n fi is the number 
density of inflatons) which together with the relation s = (2-k 2 / Ah)g*T 3 implies 
that riinfi/s = (3/4)(T r /rrii n fi). ii) Since each inflaton decays into two v c 's, 
their number density n„c = 2rii n fi which then gives n v c j s = (3/2)(T r /rrii n fi). 
iii) Supersymmetry gives an extra factor of two. 

11.2 Sphaleron Effects 

To see how the primordial lepton asymmetry partially turns into the observed 
BAU, we must first discuss the non-perturbative baryon (B-) and lepton (L-) 
number violation e3 in the Standard Model. Consider the electroweak gauge 
symmetry SU(2)l x U(1)y in the limit where the Weinberg angle 9w = and 
concentrate on SU{2)l (inclusion of 9w ^ does not alter the conclusions). 
Also, for the moment, ignore the fermions and higgs fields so as to have a pure 
SU(2)l gauge theory. This theory hasc3 infinitely many classical vacua which 
are topologically distinct and are characterized by a 'winding number' n € Z . 
In the 'temporal gauge' (Aq = 0), the remaining gauge freedom consists of time 
independent transformations and the vacuum corresponds to a pure gauge 



where g is the SU(2)l gauge coupling constant, x belongs to 3-space, i =1,2,3, 
g{x) € SU(2)l, and g(x) — > 1 as | x \— > oo. Thus, the 3-space compactifies 
to a sphere S 3 and g(x) defines a map: S 3 — > SU{2)l (with the SU(2)l 
group being topologically equivalent to S 3 ). These maps are classified into 
homotopy classes constituting the third homotopy group of S 3 , ^(S 3 ), and 
are characterized by a 'winding number' 



The corresponding vacua are denoted as | n), n G Z. 

The tunneling amplitude from the vacuum | n_) at t = — oo to the vacuum 
| n + ) at t = +oo is given by the functional integral 




(72) 



Ai = - d,g{x)g 1 {x) 

9 



(73) 




tr (d^g-^^djgi^g-^^dkg^g- 1 ^)) . (74) 




(75) 
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over all gauge field configurations satisfying the appropriate boundary condi- 
tions at t = ±00. Performing a Wick rotation, xq = t — > —ix 4 , we can go 
to Euclidean space-time. Any Euclidean field configuration with finite action 
is characterized by an integer topological number known as the Pontryagin 
number 



9 2 Im 



d x tr 



(F^F^) , (76) 



* I6tT 2 

with 2,3,4 and — ^e^vXpF p being the dual field strength. But 



\x(]?i iV F } = d^J^, where J M is the 'Chern-Simons current' given by 
= e^uap tr U v F a P - 2 -gA v A a . 



(77) 



In the 'temporal gauge' (A = 0), 



-f^r / d 4 x d»J^ = A /d 3 x J 

167T 2 7 167T^ £c 4 = ±oo J 

d 3 x e ijk tr {d l gg' 1 d j gg' 1 d k gg- 1 ) = n+ - n_ . (78) 

This means that Euclidean field configurations interpolating between the vacua 
I n+), I n_) at X4 = ±00 have Pontryagin number q = n+ — ri- and the path 
integral in Eq.[75] should be performed over all these field configurations. 
For a given q, there is a lower bound on S(A), 

S(A) >%\q\ , (79) 

which is saturated if and only if F^ u = ii 7 ^,, , i.c, if the configuration is self-dual 
or self-antidual. For q=l, the self-dual classical solution is called instanton e3 
and is given by (in the 'singular' gauge) 

A (x) - ^{x-zf m) 

AaAX) ~ g(x-zf {x-zf + f ' (80) 

where rj aflL , (a=I,2,3; /i,f=l,2,3,4) are the t' Hooft symbols with r) ai j — e ai j 
(ij'=l,2,3), r/ a4i = ~S ai , r\ al4 = S al and i] a4:4 = 0. The instanton depends 
on four Euclidean coordinates (its position) and its scale (or size) p. Two 
successive vacua | n), | n + 1) are separated by a potential barrier of height 
oc p^ 1 . The Euclidean action of the interpolating instanton is always equal to 
8-7T 2 /g 2 , but the height of the barrier can be made arbitrarily small since the 
size p of the instanton can be taken arbitrarily large. 
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We now reintroduce the fermions into the theory and observe cB that the 
B- and L- number currents carry anomalies, i.e., 

« = « = -%yf^ . (81) 

where n g is the number of generations. It is then obvious that the tunneling 
from | n_) to | rt+) is accompanied by a change of the B- and L- numbers, 
AB = AL = -n g q = -n g (n + - n_). Note that i) A(B — L) = 0, and ii) 
for q=l, AB = AL = — 3 which means that we have the annihilation of one 
lepton per family and one quark per family and color (12-point function). 

We, finally, reintroduce the Weinberg-Salam higgs doublet <j> with its vev 
given by 

<cf)>=^= ^Jj , 246 GeV. (82) 

It is then easy to see that the instanton ceases topExist as an exact solution. 
It is replaced by the so called 'restricted instanton'C3 which is an approximate 
solution for p <C v^ 1 . For | x — z |<C p, the gauge field configuration of the 
'restricted instanton' essentially coincides with that of the instanton and the 
higgs field is 

v ( {x-zf \ 1/2 fO 



[J^wf?) [i) ■ < 83 > 

For | x — z |^> p, the gauge and higgs fields decay to a pure gauge and 
the vev in Eq.^2] respectively. The action of the 'restricted instanton' is 
S r i = (87r 2 /g 2 ) + ir 2 v 2 p 2 + ■ ■ •, which implies that the contribution of big size 
'restricted instantons' to the path integral in Eq.f75| is suppressed. This justi- 
fies a posteriori the fact that we restricted ourselves to approximate instanton 
solutions with p <C v^ 1 . 

The height of the potential barrier between the vacua | n), | n+ 1) cannot 
be now arbitrarily small. This can be understood by observing that the static 
energy of the 'restricted instanton' at x^ — (A is the higgs self-coupling), 

E b (p)K^\ + L*v 2 p 2 + ^ 2 vV, (84) 



9 P 



is minimized for 




(85) 
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and, thus, the minimal height of the potential barrier turns out to be E n 
Mw/a-w ( Mw is the weak scale and aw — g 2 /^)- The static solution 
which corresponds to the top (saddle point) of this potential barrier is called 
sphaleronEJ and is given by 

i=t ) Mfxr^=^/ 1 (C)f-f(j) 1 (86) 

where £ = 2M]yr, f is the radial unit vector and the 3-vector f consists 
of the Pauli matrices. The functions /(C), which can be determined 

numerically, tend to zero as £ — ► and to 1 as £ — * oo. The mass of the 
sphaleron is estimated to be 

E sph = k, 1.5 < k < 2.7 , for < A < oo , (87) 

aw 

and lies between 10 and 15 TeV. 

At zero temperature the tunneling from | n) to | n + 1) is utterly sup- 
pressed l£| by the factor exp(— Sir 2 /g 2 ). At high temperatures, however, ther- 
mal fkictuations over the potential barrier are frequent and this transition can 
occurEil with an appreciable rate. For Mw ~ T S T c (T c is thejCritical temper- 
ature of the clcctroweak transition) , this rate can be calculated^ by expanding 
around the sphaleron (saddle point) solution and turns out to be 

v(T) 9 

T « 10 4 n g exp(-£ sph (T)/T) . (88) 

Assuming that the electroweak phase transition is a second order one, v(T) and 
E sph {T) oc (1-T 2 /T 2 ) 1 / 2 . One can then show that T > H for temperatures T 
between ~ 200 GeV and ~ T c . Furthermore, for temperatures above T c , where 
the sphaleron solution ceases to exist, it was arguedcd that we still have r 3> H. 
The overall conclusion is that non-perturbative B- and L- number violating 
processes are in equilibrium in the universe for T > 200 GeV. Remember that 
B — L is conserved by these processes. nr-i 

Given a primordial L- number density, one can calculate I13H the result- 
ing ub/s (tib is the B- number density). In MSSM, the SU(2)l- instantons 
produce the effective operator (in symbolic form) 

2 = (qqql) n °(hWhW)W 4 , (89) 

and the SU(3) C - instantons the operator 

3 - (qqu c d c ) na 9 6 , (90) 
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where q, I are the quark, lepton SU(2)l- doublets respectively, u c , d c the up, 
down type antiquark SU(2)l- singlets respectively, g, W the gluons and W- 
bosons and tilde represents their superpartners. We will assume that these 
interactions together with the usual MSSM interactions are in equilibrium at 
high temperatures. The equilibrium number density of ultrarelativistic parti- 
cles An = n part - n ant i par t is given by 



An 15(/ 



s 47T 2 g» 



where g is the number of internal degrees of freedom of the particle under 
consideration, /z its chemical potential and e = 2 or 1 for bosons or fermions 
respectively. Each interaction in equilibrium implies that the algebraic sum 
of the chemical potentials of the particles involved is zero. Solving these con- 
straints, we end up with only two independent chemical potentials, fi q and fig, 
and the baryon and lepton asymmetries can be expressedCJ in terms of them 
as follows ^ 

IT = 47r 2 g,r ^ 6 " gAtg ~ ~ 9 ^ ' 

where fi(n 9 ) i s a knownEl function. Now soft supersymmetry breaking cou- 
plings come in equilibrium at T < 10 7 GeV since their rate Tg i=s m 2 s /T > H w 
30 T 2 /Mp (ms is the soft supersymmetry breaking scale). J-n particular, the 
non-vanishing gaugino mass implies \ig = and Eqs.|9^ givec3 

ns = 4(1 + 2n g ) n B -L 
s 22n g + 13 s 

Equating tib-l/s with the primordial til/s, we have ns/s = (— 28 / 79) (til/ s), 
for n g = 3. 



12 Supersymmetric Hybrid Inflation and Neutrinos 

The supersymmetric model based on the left-right symmetric gauge group 
Glr, which we discussed in Sec.^, provides a framework within which hy- 
brid inflation, baryogenesis and neutrino oscillations arel£3 closely linked. This 
scheme, supplemented by a familiar ansatz for the neutrino Dirac masses and 
mixing of the two heaviest families and with the MSW resolution of the solar 
neutrino puzzle, implies EJ relatively stringent restrictions on the mass of the 
tau-neutrino, m„ T , and the /j, — t mixing angle, # MT . 
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In order to obtain information about the light neutrino masses and mixing 
of the two heaviest families, we ignore the first family assuming it has small 
mixings. The relevant 'asymptotic' (at the GUT scale Mx) 2x2 mass matrices 
are M l , the mass matrix of charged leptons (l c , I), M D , the Dirac mass matrix 
of neutrinos (y c : , v), and M R , the Majorana mass matrix of v c 's. We shall 
first diagonalize M l , M D : 

M' -» M" = U l °M l U l = ( 771,1 ] , (94) 
M D —* M D ' = U vC M D U v = ( m ° D \ (95) 

\ m 3 / 

where the diagonal entries are positive. This gives rise to the 'Dirac' mixing 
matrix U h ,j( U l in the leptonic charged currents. Using the remaining phase 
freedom, we can bring this matrix to the form 

U^U^f C0S6 ° S[ *° D D ) , (96) 
\ — sinfl cos0 / 

where 9 D (0 < 9 D < n/2) is the 'Dirac' (not the physical) mixing angle in 
the 2 — 3 leptonic sector. In this basis, the Majorana mass matrix can be 
written as M R = U^MoU' 1 , where M = diag(M 2 , M 3 ), with M 2 , M 3 (both 
positive) being the two Majorana masses, and U is a unitary matrix which can 
be parametrized as 

- ( C ° S9 S[U9 C ~ iS \ ( ^ rq _s 



-sm9e l ° cos8 

with < 9 < n/2 and < S < ir. The light neutrino mass matrix is 

- ^ 1 / P iCt2 \ / P^ 2 \ 

™ = - MD ' MH MD '={ e iasJ^^)( e <«.). 08) 

where ^(9,6) depends also on M 2 , M 3 , mf, m^. 

We will denote the two positive eigenvalues of the light neutrino mass 
matrix by m 2 (or m Vfl ), m 3 (or m„ r ). Recall that all the quantities here 
(masses, mixings) are 'asymptotic'. The determinant and the trace invariancc 
of rrvm provide us with two constraints on the (asymptotic) parameters: 



P m P\ 



m^ra 



2 



777 2 7713 = V M 2 M 3 ' (99) 
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(mg 2 c 2 + mf V) 2 2jmf 2 - 2 ) 2 cV cos 2j 
M 3 2 + M 2 M 3 

where 8, 6 are defined in Eqj97[ c = cos6>, s = sin 6. Note that the phases 
«2 1 in Eq.|98| cancel out in these constraints and, thus, remain undetermined. 
The mass matrix m is diagonalized by a unitary rotation V on v 's: 

v = ( elP2 \ ( cosip sin ^ e ~ l£ ^ noil 

^ J \ -sin^e" cos^ J ' 1 ' 

where < ip < ir/2 , < e < 7r. The 'Dirac' mixing matrix in Eq.^ is now 
multiplied by on the left and, after phase absorptions, takes the form 

cos# 2 3 sin# 23 e-^ \ 
-sin6» 23 e iS23 cos6» 23 J 1 { U > 

where < #23 < tt/2 , < 623 < 7r. Here, #23 (or MT ) is the physical 
mixing angle in the 2 — 3 leptonic sector and can be determined from the 
requirement that its cosine is equal to the modulus of the complex number 
cos^costf^ + sini J 9sin#' D e 1 ^ -6 ) , where — n <£ — e = /?2 — P3 — e<vr. The 
phases [3 2 , fiz and £ remain undetermined due to the arbitrariness of a 2 , a 3 . 
Thus, the precise value of #23 cannot be found. However, we can determine 
the range in which 6> 23 lies: | tp - 6 D \ < 6 23 < if + 9 D , for ip + 6 D < tt/2 ■ 

We now need the asymptotic values of , D . Approximate SU(4) C - 
invariance in the up quark and neutrino sectors gives m 2 = m c , = m t , 
s'm6 D =| V c b I 'asymptotically'. Renormalization, of light neutrino masses 
and mixing between Mx and Mz is also included E3E 2 ] assuming the MSSM 
spectrum and large tan/3 ss rut/rub- In the framework of 'hierarchical' light 
neutrino masses (m 3 3> to 2 wi), the small angle MSW resolution of the 
solar neutrino puzzle implies^ 1.7 x 10~ 3 eV 5 ni2 ;$ 3.5 x 10~ 3 eV. Finally, 
m 3 is restricted by the cosmologicaL bound m 3 < 23 eV (for h ss 0.5). 

We are now ready to derive oEIl useful restrictions on M2.3- Assume that 
both M2,3 < rrii n fi/2. Then the inflaton predominantly decays to the heaviest 
of the two. The determinant condition implies that the lowest possible value 
of the heaviest Af 2 , 3 is aboutJLO 11 GeV giving T r > 10 10 GeV, in conflict 



with the gravitino constraint El (see SecJlOl). So we are obliged to require 
that 1.72 x 10 13 GeV w m inf i/2 < M 3 < 2.5 x 10 13 GeV, where the upper 
bound comes from the requirement that the coupling constant of the non- 
renormalizable term responsible for the mass of the heaviest v c does not exceed 
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Figure 1: The allowed regions in the m Vfl , m„ T plane. 



unity. Th is requirement also implies that xq S 2.6, which makes our choice 
in Sec. 9/2, xq w 2, central. In summary, we see that i) M3 is constraint in a 
narrow range, and ii) the inflaton decays to the second heaviest right handed 
neutrino, v c , with mass M%. 

Baryons can be produced, in the present scheme, only via a primordial 
leptogenesis from the decay of v c 's emerging as decay products of the inflaton 
(see Sec. 11. f). The lepton asymmetry is then partially converted into the ob- 



served baryon asymmetry of the universe by 'sphaleron' effects as we explained 
From Eq.fn], the lepton asymmetry is 



in Sec.11.2 



til 
s 



9T r 



M 2 



sin 26 (m^-mf 2 ) 2 



87rm in// M 3 !(/*«) | 2 (m 



D 2 2 



2 c' 



) 



(103) 



Renormalization effects should also be included in this formula. Assuming the 
MSSM spectrum between 1 TeV and Mx, we saw in Sec. 11.2 that the observed 
baryon asymmetry ns/s is related to tll/s by ub/s = (— 28/79) (nx,/s). 

We will now extract Lj restrictions on light neutrino masses and mixing. 
Take a specific value of M3 (in practice, we take its two extreme values mj n /;/2 
or 2.5 x 10 13 GeV). For any pair 7712, m^, we use the determinant condition 
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Figure 2: The allowed regions in the v^-v T oscillation plot. 



to evaluate and, thus, T r . The gravitino constraint El (T r < 10 9 GeV) 
then gives a lower bound in thn rri2, plane. This bound together with the 
small angle MSW restriction E3 on 1112 yields a lower bound on 7713, namely 
m 3 > 0.9 eV (for M 3 = 2.5 x 10 13 GeV) or 1.3 eV (for M 3 = m m/i /2). 
The trace condition is so lved with respect to S = 6(9), < 9 < n/2, which 
is then substituted, in Eq jl03| to yield ul/s = (nz/s)(9). Imposing the 'low' 
deuterium boundEflon ns/s (0.02 < fish 2 S 0.03), we find the range of 9 where 
this bound is satisfied. If such a range exists, we keep m2, m 3 as satisfying 
the baryogenesis constraint. This gives an upper bound in the 777,2, plane 
which together with the MSW restriction on 7772 yields an upper bound on 7773, 
namely 7773 £ 5.1 eV (for M 3 = 2.5 x 10 13 GeV) oriL8 eV (for M 3 = m inf i/2). 
The allowed area in the m,2, 1TI3 plane is depicted Ej in Fig|l], where the thick 
solid (dashed) line corresponds to_/Vf 3 = m in fi/2 (M3 = 2.5 x 10 13 GeV). The 
overall allowed range for m„ T is Ej 1 eV Z m„ T ^ 9 eV, which is interesting 
since the value of m VT required by the cold plus hot dark matter scenario 
for large scale structure formation in the universe is centrally located in this 
range. 

The mixing angle 9^ T can also be restricted. For every allowed 7772 , 7773 pair 
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and every 9 satisfying the baryogenesis constraint, we construct V in Eq. 101 



and, consequently, tp, e and the allowed range of (9 pr , | ip — 9 D \< 9^ T < ip + ffl 7 . 
The union of all these ranges for all allowed 9 's and mi 's for a given gives 
the range of 9^ T which is allowed for this value of 7713 . All these ranges for all 
allowed ma 's constitute the allowed area on the oscillation diagram, which is 
depictedEJ in Fig.|| (notation as in Fig.|l|) in confrontation to past, ongoing and 
planned experiments. The central part of this allowed area will be tested by the 
ongoing Short Baseline Experiments (SBLE) at CERN, NOMAD/CHORUS. 
Possibly negative result from NOMAD/CHORUS will exclude a significant 
part of the allowed domains in Figs|l], || reducing the upper bound on the 
tau-ncutrino mass, m Ur , to 3.7 eV. The new CERN SBLE(TOSCA) together 
with the new CERN Medium Baseline Experiment (MBLE) ICARUS-JURA 
(600t,DIS) will cover all our predicted area on the oscillation diagram. 

In summary, hybrid inflation, baryogenesis and neutrino oscillations have 
been linked in the context of a supersymmetric model based on a left-right sym- 
metric gauge group. Our scheme leads to stringent restriction on txi Vt and 9^ 
to be tested by ongoing and planned experiments. These restrictions are de- 
rived by mainly 'physical' arguments (gravitino and baryogenesis constraints) 
supplemented by a 'minimal' input from fermion mass matrix ansaetze (only 
3 input parameters ) and experiments (MSW resolution of the solar neutrino 
problem). The choice of the gauge group is crucial since, in this case, </> has 
the quantum numbers of v c and , thus, decays to v c 's producing an initial lep- 
ton asymmetry. As a consequence, the gravitino and baryogenesis constraints 
restrict the neutrino parameters. 
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